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One- and two-dimensional harmonic polylogarithms, HPLs and GPLs, appear in calculations of multi-loop 
integrals. We discuss them in the context of analytical solutions for two-loop master integrals in the case of massive 
Bhabha scattering in QED. For the GPLs we discuss analytical representations, conformal transformations, and 
also their transformations corresponding to relations between master integrals in the s- and f-channel. 



1. Introduction 

Multiloop calculations may be treated with nu- 
merical or analytical approaches. The latter is 
limited to issues with few, typically at most three 
different scales. But it has well-known attractive 
features, notably the knowledge of the explicit 
structure of the singularities and a good control 
of numerical stability. 

The analytical approach to be discussed here 
relies on differential equations |ll2j . Iterative so- 
lutions in powers of e = (4 — D)/2 are inher- 
ently connected with certain classes of special 
functions. For two-scale problems (e. g. QED 
self-energies and vertices), the harmonic polylog- 
arithms H{{a},x) (HPLs) have been introduced 
[3]; the index vector {a} has elements 1, 0, —1. We 
just mention that HPLs with some generalized ar- 
guments are introduced in A generalization of 
HPLs for three-scale problems (e.g. QED boxes) 
are the two-dimensional harmonic polylogarithms 
G{{b},x) (GPLs); the index vector {b} has ele- 
ments 1,0, —1, but now also those depending on 



a second kinematic variable y. This was observed 
in 5 and systematically worked out in 6, for a 
planar massless problem. There one may find an- 
alytical expressions for the GPLs until weight 3 
with indexes 0, 1, z,l — z. Another generalization 
is performed in T . Algebraic non- linear factors 
arc introduced as being motivated by the physical 
nature of the problem treated there. 

Let us shortly review some facts on HPLs. 
HPLs up to weight 4 are expressed, with few 
exceptions, in terms of Nielsen Polylogarithms 
ITSU' . The exceptions are of weight 4 and may 
be expressed by few integrals 
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where Ii{x) — H[—l,0,0,l,x]. However, it can- 
not be excluded that some relations among them 
exist. In fact, a relation between 1^ and I4 exists 
|1U) . In fact, exploiting that J dy (y / (1 — y)) is 
known to Mathematica (after substituting y/{l — 
y) —>■ z) and is related to /a — 1^, one may derive 
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As far as numerical solutions are concerned, 
the problem is solved for any weight of HPLs. 
The numerical evaluation of HPLs with indexes 
{0,1,-1} (and the analytic continuation) is de- 
scribed in jll| and may be performed with the 
Fortran program hplog. These HPLs have been 
programmed also in Mathematica [Q. The sys- 
tematic numerical treatment of GPLs with in- 
dexes {0, 1, 1 — z, —z} is given in and per- 
formed with Fortran program tdhpl.^ Please no- 
tice a change of notation compared to (H] as de- 
scribed in the Appendix of |13| . 

In this paper, we will discuss features of the 
GPLs which we are exploiting for a study of mas- 
sive QED two-loop boxes. These GPLs were used 
in e.g. |14ll5j and |9llfill7| . In the next section 
we describe the general idea of solving differen- 

^This is in fact not sufficient to cover our physical cases; 
we need indexes {0, 1, —1, 1 — z, —z}). 



tial equations for Feynman integrals with GPLs, 
then we sketch a simple one loop case, and what 
finally follows is a section on some algebraic rela- 
tions connected with GPLs as we are using them 
for Bhabha scattering. 

2. Differential equations 

Master integrals (Mis) M may be determined 
as solutions of appropriate differential equations. 
The general idea is simple.^ Let M fulfill some 
differential equation with respect to x. 



dx 



M{x) = A{x)M{x) + B{x), 



(6) 



where x is usually the external scale like s (or t) 
or a conformal counterpart x (or y); for defini- 
tions and details see The inhomogeneity B is 
a linear combination of Mis of lower complexity 
and assumed to be known from an earlier stage 
of iteration. For simplicity we assume here M to 
depend only on x. If H{x) is a solution of the 
homogeneous equation. 



dx 



H{x) = A{x)H{x), 



then the full solution is given by 



M{x) = H{x) 
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The solutions in arbitrary dimension D are gener- 
ally combinations of generalized hypergeometric 
functions which are difficult (or at least tedious) 
to find and to expand in powers of e. This state- 
ment is true already for the massive QED one- 
loop box 19' , where Appell hypergeometric func- 
tions and a Kampe de Feriet function appear. An 
alternative idea is to expand © in e: 
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where a is fixed by the physical problem and /3 
chosen reasonably. Then one may try to solve 



nice pedagogical introduction is 1181 . 
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iteratively the system of equations for the (x) : 



Q — Z 



(11) 



—m\x) = 2_^a^{x) TO* \x) + b''{x). 

Let us mention that A should not be singular in 
e (i > in (|10|l 'l. Otherwise some m*(x) on the 
LHS of Hll|l would depend on a higher order com- 
ponent TO*~^ (x) of M sitting on the RHS, and the 
recursion could not be solved. With this assump- 
tion the solution is of the form 



to'(x) = H{x)i Const + 
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H{x) is the solution of the homogeneous equation 
( |7| l and is the same for all orders in e. 

In general, some sets of Mis will fulfill a sys- 
tem of linear differential equations. For massive 
Bhabha scattering so far only one such system 
B514m, the two-loop boxes with 5 lines, 4 of them 
massive, with two Mis has been fully solved ana- 
lytically in the language of HP Is and GPLs |9I15| . 
Another analytically fully solved 4-point MI is 
B512ml 16 . We will comment on the more com- 
plex system of 4 Mis B512m3 in the context of 
GPLs in Section H 

3. The Bhabha one-loop box: An illustra- 
tion 

A differential equation for the one-loop box 
B412m is 3 
d 

s— B412iii[s, t] = a B412m[s, t] + b SE212m[s] 

OS 

+c Tlllm + d SE210m[t] + e V311m[s]. (13) 

The tadpole T, self energies SE, and vertex V are 
known. The factor a of the homogeneous part is 

_ 8 + s'^ -2t + s{-6 + t + et) 
" ~ {-4 + s){-A + s + t) ' 
After a change to conformal variables 
^-s{t)+A- V-s(i) 



x(y) = , , 

^-s{t)+4 +^/^' 

^For notations see 



(14) 



the a has a simple rational denominator: 

1 1 1 „ r 1 

" (l + x)^ (x + y) (l + xy) ^ 

Such an expression can now be decomposed into 
terms with monomial denominators depending on 
X only. 
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and those which depend both on x and y: 



(16) 



X + y I + xy 

Factors b, c, d, e in 1)13(1 follow the same structure. 

The monomials become the kernels for two spe- 
cific classes of GPLs: 

G{b,{a},x) = f dx'g{b,x')G{{a},x'), (17) 
Jq 

b = -y,-l/y. 

Analogous relations define the HPLs, now with 
the index vector {0,1,-1}. If a simpler master 
was expressed by a structure like (|17|) . then m* 
will be easily determined by iterating the GPLs. 

In this way, after performing the e expansion of 
the objects in p3(l . the MI B412m can be solved 
systematically step by step to the desired order 
in e. For more details on the specific example we 
refer to [15j . 

4. Algebra of GPLs 

4.1. Analytical representation of GPLs 

Up to weight two, it makes not much effort to 
find analytical representations for GPLs in terms 
of Nielsen polylogarithms. If we consider only 
the index vectors {—y, — l/j/} and {0, —1, 1}, we 
have e.g. at weight two 25 GPLs (see Table 
of which we can choose freely 5 as the irreducible 
integrals. 

We put some results on GPLs up to weight four 
in file GPL . m at [20] . With conventions we follow 
|15ll6ll7l^ . Starting from weight three some of 
the GPLs in GPL . m are written as numerical inte- 
grals, e.g. 



G[-y, -y, -1,2;] 



G[-y-l,z]/iy+z)dz.{18) 



4 



Similarly to the case of HPLs discussed in the In- 
troduction, some of these integrals can be surely 
solved further to their analytical form. 



Table 1 

Number of GPLs for the index vectors 
{— 2/, — 1/y} and {0,-1,1}. The second row 
counts HPLs with indexes {0, —1, 1}. 



weight 1 


weight 2 


weight 3 


weight 4 


2 GPLs 

3 HPLs 


16 GPLs 
9 HPLs 


98 GPLs 
27 HPLs 


544 GPLs 
81 HPLs 



4.2. Interchange of arguments in GPLs 

Sometimes in the course of solving Mis we want 
to use the knowledge of differential equations in 
a fixed channel, but in both s- and i-operators 
(or equivalently in x and y). As an example may 
serve the MI B512m3d3. The homogeneous part of 
its equation, if derived with the s-operator, van- 
ishes, so that its solution is a constant. Then, the 
solution of the inhomogeneous equation is easy 
but it is not possible to find the constant term 
(a function Cy{y)) in the usual way by exploiting 
the knowledge of analyticity. However, we can 
alternatively determine the same function from 
an equation derived with the i-operator, having 
now another constant term (a function Cx{x)). 
An easy way to determine Cy(i/) is to compare 
the functional dependences of both solutions, e.g. 
in the limit x = I, and to end up with an un- 
known true constant Cx{^)- It depends on poly- 
nomials in C2,C3j • ■ • The Cx{i), in turn, can be 
numerically fitted by knowledge of B512m3d3 at 
some Euclidean kinematic point using the sec- 
tor decomposition method |21l22j . In this way, 
B512m3d3 has been determined up to e" (yet un- 
published). Sometimes it may appear from the 
very beginning more convenient to solve some dif- 
ferential equations in the t-operator (where e.g. 
functions G[—x, y] appear) and rewrite this so- 
lution later for the s-operator where other Mis are 
to be solved (and where functions G[—y, x] ap- 
pear). The interchange of arguments can also be 



used to find limits of GPLs when e.g. a; = 1 or 

Useful relations for the case of massive Bhabha 
scattering are tabulated up to weight four in 
the file GPLtransf .m at HOI- AU of them are 
based on the type of identities described in 
For instance, the first item of weight two in 
GPLtransf .m, 

Ghy,-l,x] = H[-l,x]H[l,y] 
- H[0,x]H[l,y] 
+ H[-l,x]H[~x,y] 
+ H[0,-l,x] 

+ H[l,0,y]-G[l,-x,y] (19) 
can be solved using the relation 
H[m{z);y] = H['m{z = 0);y] 

+ ^ dz'^H[m{z');y] (20) 

after interchanging the differential and integra- 
tion operations, combined with a decomposition 
into basic polynomials. 

4.3. Conformal transformations 

There are two different definitions in use for 
conformal transformations x (y) of s (t) in mas- 
sive QED. These are (|14|l . being used e.g. in 
prci . and 



x'iy' 



1 



(21) 



being used e.g. in j23l24| . However, x and x' are 
connected by a relation which is also used for a 
transformation of variables in HPLs 0: 



1 



y 



l + y' 



(22) 



So, if we have to transform Mis defined in terms 
of X to Mis defined in terms of x', we have to know 
the corresponding relations between GPLs.^ Here 
is an example for this: 



G{-y,0,x) 



H[-y, 0,1] 



dz 



H[0, z] 

JEZ 

i+y' 



z + 



* Results given in I23I24I are written directly using Nielsen 
Polylogarithms. 
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[ dz[H{l,z) + H{-l,z)] 
Jo 



1 



(23) 



1 + z 1 — zy' 
Interchanging arguments, 

H[-y,0,l] = -C2 + i^[0,-l,y]- i/[0,0,2/], 

using conformal transformations for HPLs 
and integrating the expression in (|23|) we finaUy 
get 

G[-y,0,x] = ~C2-H[-l,l,y']-H[l,l,y'] 
^ H[-l,l]iH[-l,y']+H[l,y']) 

+ i?[0, -l,l]-iJ[0,0,l] 

+ H[-l,-l,x']+H[-l,l,x'] 

+ G[l/y',-l,x']+G[l/2/',l,x']. 

(24) 

Let us note that the conformal transformation 
(j22|l extends the original set of GPLs by those 
with additional arguments {+y,+l/y}. Again, 
results are given in a file GPLconf .m at |2()j . 

In conclusion, we have discussed some basic 
features of GPLs which are used in an ongoing 
determination of fully analytical results in mas- 
sive Bhabha scattering in QED. For other inter- 
esting properties of GPLs we have to refer to the 
literature; e.g. analytical continuation of GPLs 
is explored in JT]. A complete GPL package for 
general use in Mathematica will be given else- 
where. 
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